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ABSTRACT 

Simple pure luminosity evolution (PLE) models, in which galaxies brighten at high 
redshift due to increased star- formation rates (SFRs), are known to provide a good 
fit to the colours and number counts of galaxies throughout the optical and near- 
infrared. We show that optically defined PLE models, where dust reradiates absorbed 
optical light into infrared spectra composed of local galaxy templates, fit galaxy counts 
and colours out to 8|xm and to at least z ~ 2.5. At 24-70 |xm, the model is able to 
reproduce the observed source counts with reasonable success if 16% of spiral galaxies 
show an excess in mid-IR flux due to a warmer dust component and a higher SFR, in 
line with observations of local starburst galaxies. There remains an under-prediction of 
the number of faint-flux, high-z sources at 24|j.m, so we explore how the evolution may 
be altered to correct this. At 160|J.m and longer wavelengths, the model fails, with our 
model of normal galaxies accounting for only a few percent of sources in these bands. 
However, we show that a PLE model of obscured AGN, which we have previously 
shown to give a good fit to observations at 850|i.m, also provides a reasonable fit to 
the -ffersc/ie//BLAST number counts and redshift distributions at 250-500|i.m. In the 
context of a ACDM cosmology, an AGN contribution at 250-870|i.m would remove 
the need to invoke a top-heavy IMF for high-redshift starburst galaxies, although the 
excellent fit of the galaxy PLE model at shorter wavelengths would still need to be 
explained. 

Key words: infrared: galaxies - dust, extinction - galaxies: evolution - submillimetre 



1 INTRODUCTION 

The process by which galaxies formed and evolved into the 
systems we see today remains one of the most important 
questions in modern astronomy. We do not currently have a 
comprehensive picture of how and when most of the stellar 
content of the universe was assembled. 

In studies of galaxy evolution there are, broadly, two 
ways to proceed. One is the ah initio approach, in which 
galaxy evolution is modelled using a chosen prescription for 
processes including the collapse and merging of dark mat- 
ter halos, gas heating and cooling, star-formation and feed- 
back (e.g. Cole et al., 2000). The other approach is more 
observationally led, beginning with a description of the lo- 
cal luminosity functions (LFs) and spectral energy distri- 
butions (SEDs) of galaxies and 'winding the clock back' to 
build up a picture of the high-redshift universe (Bruzual 
& Kron, 1980; Koo, 1981; Shanks et al., 1984; Pozzetti et 
al., 1998). These approaches are often called 'forward evo- 
lution' and 'backward evolution' respectively (for a review, 
see §5.2 of Hauser & Dwek, 2001). While forward evolution 
models can potentially offer greater insight into the physical 
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processes that may be driving galaxy evolution, the advan- 
tage of backward evolution models is that they begin with 
a phenomenological description of the local universe, which 
is comparatively well understood and well constrained. 

The semi-analytical models used to simulate structure 
formation in the ACDM cosmology employ a forward evo- 
lution, ah initio approach. In these models, galaxy forma- 
tion is a hierarchical process, in which stellar mass is built 
up over extremely long timescales through the merging of 
smaller structures (Cole et al., 2000). However, hierarchical 
models have faced significant observational challenges. To 
summarise a wide range of studies, the broad observational 
picture is that the most massive galaxies host old stellar 
populations (Franx et al., 2003; Cimatti et al., 2004; Glaze- 
brook et al., 2004; Thomas et al., 2005; Bernardi et al., 2006) 
apparently indicating a high redshift of formation, and show 
little evolution in their number density since z « 1 (Wake et 
al., 2006; Brown et al., 2007; Cool et al., 2008), suggesting 
no significant dynamical evolution in the last several Gyr. 
These observations are contrary to the expectation from hi- 
erarchical evolution, where massive structures would be ex- 
pected to have assembled since 2; « 1. 

The observations are more in line with the predictions of 
passive evolution models, which assume that galaxies formed 
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early and were in place at some high redshift, thereafter 
evolving to the present day only in terms of a decreasing 
star-formation rate. Instead of a hierarchical, merger-driven 
process of structure formation, passive evolution models as- 
sume galaxies formed from massive gas clouds in compara- 
tively short-lived 'monolithic collapse' events. 

In this paper we use a model which employs the back- 
ward evolution approach and pure luminosity evolution, be- 
ginning with a local LF and evolving it to higher redshift us- 
ing k- and/or (k -\- e)-corrections with no density evolution. 
We model galaxy evolution by simply assuming that star- 
formation declines at different rates in early- and late-type 
galaxies. This model has been used extensively in the past 
at optical and near-infrared (NIR) wavelengths (Metcalfe et 
al., 1995, 2001, 2006; McCracken et al., 2000) where, with 
modifications (see Metcalfe et al., 2006), it has been success- 
ful in matching observations from the U to the K band and 
out to z « 4. We now test the model's predictions against 
observations in the mid-IR, far-IR and sub-mm. 

These long wavelength regimes are important to galaxy 
evolution studies. Observations with the COBE satellite in 
the 1990s confirmed the existence of a cosmic infrared back- 
ground (CIB) which peaks in the far-IR at A ~ 100-300[J.m 
(Puget et al., 1998; Fixsen et al., 1998; Hauser et al., 1998). 
It is now clear that the CIB, which arises from the re- 
emission by dust of absorbed radiation, has a total intensity 
roughly equal to that of the extragalactic background at 
optical wavelengths. Since infrared sources contribute only 
~ 30% of the total energy output of the local universe (Soifer 
& Neugebauer, 1991), there is expected to be significant evo- 
lution in the populations which contribute to the infrared 
background. It is clear, then, that an understanding of the 
evolution of galaxies at long infrared wavelengths is impor- 
tant to studies of the star-formation and, since AGN may 
contribute to this background, accretion history of the uni- 
verse. 

Other authors have presented models which give a 
good fit to mid-IR, far-IR and sub-mm observations, using 
both hierarchical forward evolution models (e.g. Granato et 
al., 2000; Lacey et al., 2008) and phenomenological back- 
ward evolution models (e.g. Chary & Elbaz, 2001; Rowan- 
Robinson, 2001, 2009; Lagache, Dole & Puget, 2003; Lagache 
et al., 2004). 

Typically, these models are designed specifically to 
match infrared observations. The hierarchical ACDM model 
requires a top-heavy IMF to fit the FIR/sub-mm data 
(Baugh et al., 2005; Lacey et al., 2008), and the phenomeno- 
logical models tend to employ a range of different spectral 
types, including luminous starbursts, or allow for strong evo- 
lution in the luminosity functions. Here our approach dif- 
fers: all aspects of our model, including the modest amount 
of dust absorption, were defined in the optical bands. The 
surprise has been that this simple model continued to rep- 
resent the data well as it was applied at high redshifts and 
to longer, near-IR wavelengths. We now extend this much 
further, out to the mid-IR, far-IR and sub-mm regimes, fol- 
lowing Busswell & Shanks (2001). Determining where and 
how the model breaks down is instructive, since this may in- 
dicate a requirement for dynamical evolution or for another 
population to dominate the counts, e.g. AGN (see Hill & 
Shanks, 2010). 

This paper is structured as follows. In ^ we discuss 
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Figure 1. A BRI colour-colour diagram showing the predicted 
redshift tracks for the r = 9 spirals (soHd blue line) and the t = 1 
ellipticals (dashed red line). The data shown are R < 25 galaxies 
from deep Subaru imaging of the GOODS-N field. Our models 
match well the looping features which are clearly visible in the 
data. 



the datasets used; in ^ we discuss our galaxy models. In 
m we investigate galaxy colours. ^ and ij6] deal with pre- 
dictions of our PLE models for number counts at mid-IR 
and far-IR/sub-mm wavelengths respectively. In ij7] and 31] 
we present a discussion and a summary of our findings. 

Magnitudes throughout the paper are quoted in the 
Vega system unless otherwise stated. For consistency with 
the work of Metcalfe et al. we use a cosmology with Ho = 50 
km s~^ Mpc~\ rim = 0.3 and Qa = 0.7. 



2 DATA 

To test our model predictions of colours at Spitzer IRAC 
wavelengths, we make use of data from the « 200 arcmin^ 
COODS-North survey. This field is well suited since it has 
some of the deepest NIR and IRAC imaging available over 
a relatively wide area and has also been targetted in sev- 
eral spectroscopic follow-up surveys. For photometry we use 
the catalogue recently presented by Wang et al. (2010), 
who combine public Spitzer IRAC data (Dickinson et al., 
in prep.) with a new ultradeep A'-band image, using the K 
positions as priors and detecting robust IRAC counterparts. 
We get redshifts for these sources by matching to redshift 
catalogues from Cohen et al. (2000) , Cowie et al. (2004) and 
Wirth et al. (2004). We supplement this with the catalogue 
presented by Reddy et al. (2006), which includes spectro- 
scopic redshifts and A'-band and IRAC fiuxes for 388 high- 2: 
galaxies in GOODS-N. Our total redshift sample includes 
1,617 galaxies. 

Numbers counts of galaxies in the Spitzer IRAC (3.6, 
4.5, 5.8 and 8.0^m) and MIPS (24, 70 and 160M.m) bands 
have been measured across a wide flux range by surveys 
which combine observations from several flelds of differing 
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width and depth in order to achieve good sampling of both 
the faint and bright end of the source counts. At IRAC wave- 
lengths, we use the published counts of Fazio et al. (2004), 
who provide counts from three regions ranging from an ul- 
tradeep 5' x 10' field to the wide 3° x 3° Bootes field. A 
similar selection of fields was used by Papovich et al. (2004) 
and Dole et al. (2004) in the MIPS bands, with the combi- 
nation of deep and wide fields allowing them to determine 
number counts across 3 decades of fiux. A more recent survey 
by Bethermin et al. (2010) took this approach further, com- 
bining MIPS observations across 9 fields of various depths 
covering ~ 54 deg^ in total. We use all of these datasets, 
as well as our own counts from public GOODS-N imaging 
in the IRAC bands (Dickinson et al., in prep.) and at 24ji-m 
(Chary et al., in prep.) using SExtractor (Bertin & Arnouts, 
1996) on the image frames. 

We also compare our model to observations at 60jim, 
where counts are available from wide-field IRAS surveys 
(Lonsdale et al., 1990; Rowan-Robinson et al., 1991; Gre- 
gorich et al., 1995), although these observations do not reach 
the same depth achieved by the more recent Spitzer surveys. 

Finally, we use recent galaxy surveys in the sub-mm 
bands at 250, 350 and 500jim. These are the wavebands 
probed by the SPIRE detector on Herschel (Pilbratt et al., 
2010) and by Herschel's balloon-borne forerunner, BLAST 
(Patanchon et al., 2009). We make use of BLAST number 
counts (Patanchon et al., 2009) and redshift distributions 
(Chapin et al., 2010; see also Dunlop et al., 2009) from ob- 
servations covering f» 10 deg^. We use SPIRE number counts 
from the HerMES survey, currently covering ~ 20 deg^, both 
as published by Oliver et al. (2010) and also as presented by 
Dunlop et al. (2010). 



3 MODEL 

Here we give a brief summary of the galaxy evolution model 
we employ, and we review some previous results from this 
model at optical/NIR wavelengths. For a more detailed de- 
scription of the model we refer the reader to Metcalfe et al. 
(2006) and references therein. 



3.1 Galaxy types 

Our basic galaxy model assumes only two populations, 
broadly corresponding to early-types (old, red, passive sys- 
tems) and late-types (young, blue, star-forming spirals). 
Both types are modelled with exponentially decreasing star 
formation rates (SFRs), i.e. SFR(t) oc e~*''^. The early- and 
late-type galaxies are distinguished by different character- 
istic e-folding times, r: for early-types we use r = 1 Gyr, 
for late-types t — 9 Gyr. Based on this formula we use 
Bruzual & Chariot (2003; hereafter BC03) models to pro- 
duce template spectral energy distributions (SEDs) for these 
two galaxy types, assuming a Chabrier IMF. 

Assuming a bimodal galaxy population is clearly a 
simplification, but it is not without merit. Optical colour- 
magnitude and colour-colour diagrams show evidence of 
strong bimodality, with a so-called 'red sequence' and 'blue 
cloud', the former consisting primarily of bulge-dominated 
galaxies with old stellar populations and the latter of 
young, disk-dominated, star-forming galaxies (Strateva et 
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Figure 2. Infrared SEDs for three galaxy templates. Our stan- 
dard, bimodal models comprises the early-type (red dotted line) 
and late- type (blue dashed line) galaxies; the early- types are dust- 
less so we see only the long- wavelength tail of the stellar emission, 
while the late-types include dust with a schematic representation 
of the PAH features (see text). Also shown is an SED based on 
the starburst galaxy Mrk33 (green solid line), which we intro- 
duce at A ^ 24|xm in place of a fairly small fraction of the normal 
late-types. The PAH spectrum for Mrk33 comes from IRS obser- 
vations and the green datapoints show broadband photometry for 
Mrk33. All models and data are normalised to be equal at 2|j.m; 
the y-axis is in arbitrary units of XFx- 



al., 2001). Such plots can often be seen to show two curving 
tracks, and these are well matched by models using a r = 1 
and a T = 9 population (Metcalfe et al., 2001); this is clear 
from Fig. [T] 



3.2 Dust spectrum 

In order to apply this model at infrared wavelengths, a pre- 
scription for dust absorption and re-radiation is required. We 
invoke dust only for the late-type galaxies. Following Met- 
calfe et al. (2001) and Busswell & Shanks (2001), we apply a 
1/A extinction law normalised at 4500A in the B band with 
some normalisation Ab- That is, the flux at a given wave- 
length after absorption, F\, relates to the unabsorbed flux 
-^^,0 according to equation [T] 



Fx 



10" 



'Fx 



(1) 



The absorbed radiation is re-emitted by the dust in the 
infrared. The primary component of this emission is a mod- 
ified Planck function, which has some given dust tempera- 
ture and opacity. We apply an opacity law with the standard 
form K oc A~ , where a value of /? = would correspond to 
a perfect blackbody. In our model we set /3 — 1.5 (Dunne & 
Eales, 2001). 

Emission from polycyclic aromatic hydrocarbons 
(PAHs) is an important spectral feature between 3 and 13 
microns (Leger & Puget, 1984; AUamandola et al., 1985). In 
our spiral galaxy model we include a fairly schematic repre- 
sentation of the PAH features, designed to broadly resem- 



© 2010 RAS, MNRAS OOO.mflil 



M. D. mil & T. Shanks 



10° 







+ McCracken et aL (2000) 






yi Gardner etaL (1993) 






McLeod et al. (1995) 




10^ 


_ A Saracco et al. (1997) 




D Metcalfe et al. (2006) 






Vandame el al. (2001) 






Huang el al. (2001) 




10^ 


,. , Kummel & Wagner (2001) 


T* 


_ A Maihara et al. (2001) 


OJJ 


Vaisanen et al. (2001) 


B 




E Kong etal. (2006) 




• Imai et al. (2007) 


<n 




T Frith 2MASS(priv.coinm. 


o 


10^ 


_ ♦ lovino et aL (2005) 


'on 






« 






-a 


10^ 




^^ 


- 


U 






•a 






7; 






T3 


lU 


i 



10" 



10' 




10 



15 
K 



20 



25 




Figure 3. Differential number counts in the K band. The solid 
line is the prediction of the model used in this work; the dashed 
line shows the model of Metcalfe et al. (2006) who use an a; = 3 
IMF for early-types as described in the main text. 



Figure 4. The redshift distribution of sources in the K band. The 
histogram shows observational data from the K20 survey (Cimatti 
et al., 2002) down to a magnitude limit of ii" = 20. The solid and 
dashed lines arc as in Fig. |3] 



ble the spectra of galaxies observed in the Spitzer Infrared 
Nearby Galaxies Survey (SINGS; Kennicutt et al., 2003), 
presented by Draine et al. (2007). We show the IR spectra 
of our early- and late-types in Fig. [2] While the modelling 
of the PAH features in the late-types is somewhat rudimen- 
tary, we find it is sufficient to reproduce the observed colours 
and sources counts of galaxies which sample this part of the 
spectrum, so we consider it sufficient for the purposes of this 
paper. 

Fig. [2] also includes a third SED, which represents a 
starburst component modelled on Markarian 33 (Mrk33). 
This is an additional spectral type which we introduce when 
dealing with observations at A ^ 24nm (see i )5.21 where the 
motivation for including this component is discussed) . 

We assume a dust temperature of 30K, which is found to 
be a typical temperature of interstellar dust in star-forming 
galaxies (Farrah et al., 2003; Pope et al., 2006; Coppin et al., 
2008; Elbaz et al., 2010). We take Ab = 0.3 magnitudes for 
the normalisation of the 1/A dust absorption law; this is the 
value determined by Metcalfe et al. (2001) when using this 
model at optical wavelengths, and is a fairly conservative 
amount of extinction. The total integrated flux of the dust 
emission, in the form of both the PAH features and the 
blackbody, is normalised to be equal to the total absorbed 
flux. 



3.3 Initial mass function and luminosity function 

For each galaxy type, we use a luminosity function (LF) de- 
flned locally in the B band and correct it into the appropriate 
IR band using a rest-frame z = colour. The resulting LF 
parameters for the K band are given in Table [T] Then using 
k + e corrections determined by the BC03 evolution code, 
the LF is evolved back from the present day. The model is 
a simpliflcation of the model of Metcalfe et al. (2001, 2006) 



Table 1. Luminosity function parameters. These are the same 
as those used by Metcalfe et al. (2006). When shifting the LFs 
to another band, we use 2 colours, from a t = 1 Gyr model for 
early-types (E/SO and Sab) and r = 9 Gyr model for late-types 
(Sbc, Scd and Sdm) to get new values fof M* for a given band, 
as shown in Table [2] 



Type 


#7Af pc-3 


a 


M*K 


E/SO 


9.27 X 10-'' 


-0.7 


-24.92 


Sab 


4.63 X 10-4 


-0.7 


-24.78 


Sbc 


6.20 X 10--* 


-1.1 


-24.83 


Scd 


2.73 X 10--* 


-1.5 


-24.34 


Sdm 


1.36 X 10-4 


-1.5 


-23.71 



since rather than 5 independent colours, one for each type, 
we use only two, one for all of the early-types (E/SO and Sab) 
and one for all of the late-types (Sbc, Scd and Sdm). This 
is justified on the basis that the band-band colour differ- 
ences in the NIR are smaller than in the optical. The z = Q 
colours are taken from the BC03 model predictions. Table 
[2] shows the values of M* after shifting into MIR/FIR/sub- 
mm bands. On this basis, we determine predictions for the 
observed galaxy number counts and colours. 

Metcalfe et al. (2001, 2006) used a Salpeter IMF 
(Salpeter, 1955; x — 1.35) for the k + e corrections of spirals 
but used an a; = 3 IMF, cut at M > 0.5Mq, for the fc -f- e 
corrections of early-types. This was to reduce the passive 
evolution in the K-band to almost the level of the K-band 
fc-correction to prevent too many z > 1 galaxies being pre- 
dicted at K < 20, which would be the case for the Salpeter 
IMF. Here we are using the Chabrier IMF (Chabier, 2003; 
2; = 1.3 at M > IMq) with an SFR e-folding time of r = 9 
Gyr for the spirals; this produces k + e corrections very simi- 
lar to the Salpeter with r = 9 Gyr. For the early-type galax- 
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Table 2. Characteristic magnitudes, M*, for the luminosity function of our galaxy types in all of the mid-IR, far-IR and sub-mm 
wavebands used in this paper. Since the ellipticals are dustless in our model their contribution falls essentially to zero in the FIR/sub- 
mm, therefore these colours are omitted. As described later in the paper, at 2A\ua and above we replace the Scd population with a 
population based on the starburst galaxy Mrk33. Note that for consistency with Metcalfe et al. we use Vega colours out to 24^m. At 
60|a.m and above this is infeasible, so the magnitudes given for these bands are in the AB system. 



Type 



KK 



M* 



MS 



M* 



Mt 



Mi 



E/SO 


-25.08 


-25.10 


-25.17 


-25.27 


-25.36 


- 


- 


- 


- 


- 


- 


- 


Sab 


-24.94 


-24.96 


-25.03 


-25.13 


-25.22 


- 


- 


- 


- 


- 


- 


- 


Sbc 


-25.29 


-25.27 


-26.49 


-27.81 


-30.19 


-26.07 


-26.41 


-26.61 


-25.61 


-24.68 


-23.67 


-22.56 


Scd 


-24.80 


-24.78 


-26.00 


-27.32 


- 


- 


- 


- 


- 


- 


- 


- 


Starburst 


- 


- 


- 


- 


-32.92 


-27.58 


-27.74 


-27.60 


-26.69 


-25.87 


-25.00 


-23.44 


Sdm 


-24.17 


-24.15 


-25.37 


-26.69 


-29.07 


-24.95 


-25.29 


-25.49 


-24.49 


-23.56 


-22.55 


-21.44 
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Figure 5. Colour-redshift plots showing how our models compare to X < 22.5 sources from GOODS-N spectroscopic surveys. Red and 
blue lines are PLE redshift tracks for our early-type and late-type models respectively. Black circles indicate the median of the data in 
Az = 0.2 bins. Overall there is a good agreement between the data and the model predictions. 



ies, no I = 3 (or Scale, 1986) IMF was easily available from 
BC03 and so here for predicting number counts and 71(2) 
in the K band and redward, we have simply assumed the 
fc-correction from the Chabrier IMF, with r = 1 Gyr rather 
than T — 2.5 Gyr to account for the degeneracy between 
IMF slope and SFR in the k + e predictions (Metcalfe et 
al., 2001). We have checked that these are good approxima- 
tions to the a; = 3 predictions in the near-/mid-IR bands. 
Ultimately, however, the choice of IMF and SFR for early- 
types has little impact on our results since the early-types 
are taken to be dustless, and their contribution is therefore 



negligible in most of the mid-IR, far-IR and sub-mm wave- 
bands which are the focus of this paper. 



In Fig. Owe show the model predictions for the K band 
number counts generated in this way and we note the model 
fits the counts very well across the full range of magnitudes. 
We also see that the model agrees well with the 3; = 3 
model of Metcalfe et al. In Fig. |4]we show the K band red- 
shift distribution compared to the data from the K20 survey 
(Cimatti et al., 2002) where again there is excellent agree- 
ment between the model and the data, with the fit being 
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Figure 6. Colour-colour plots showing how our models compare 
to observational data from GOODS-N. The contours show _ft' < 21 
galaxies, with the lowest contour plotted where there are 2 ob- 
jects in a 0.1 X 0.1 mag bin - below this level galaxies are shown 
as individual crosses. The red (dashed) and blue (solid) lines are 
PLE redshift tracks for our early-type and late-type models re- 
spectively. The square marker on each track indicates z = 0.5 
and the circles further along the lines mark 2 = 1, 2 and 3. For 
reference, in panel (c) we also show the selection criteria for the 
AGN wedge (dashed line; Stern et al., 2005). 
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4 GALAXY COLOURS 



3 model of Metcalfe 



In Fig. [5] we present our model predictions for colour-redshift 
tracks in near-/mid-IR bands. The predictions are compared 
to observed colours of _?S' < 22.5 galaxies in the GOODS-N 
catalogue of Wang et al. (2010), described in ij2l Overall, 
the agreement between the data and our model in Fig. [5] is 
very good, suggesting that a simple optically-defined PLE 
galaxy model can match colours in the mid-IR out to at 
least A = 8ii,m and z = 2.5. The predicted spiral galaxy 
tracks in Fig. [5] show a 'bump' in K — [3.6] at z ~ 0.1 
and in A' — [4.5] at z ~ 0.5. These arise from the PAH 
emission feature at A = 3.3(i.m passing through the 3.6 and 
4.5 micron bands, respectively. In the K — [4.5] : z plot there 
is clear evidence of this bump in the data, well matched 
by the model, however the corresponding bump is not as 
clearly seen in the K — [3.6] : z plot. This PAH line is also 
responsible for the significant peak seen in the [3.6] — [4.5] 
colour at 2: ~ 0.4, which again is replicated by the data. 

The difference between the early- and late-type models 
is most pronounced in the [5.8] — [8.0] colour (Fig. [5] lower- 
right panel), because the 5.8vi.m and 8.0(i.m bands sample the 
region of the spectrum where the contribution from dust is 
important relative to that from starlight. Since our early- 
types are modelled without dust it is the predicted colours 
of the late-types which give a better fit to most of the data 
here. 

In Fig.|6]we show colour-colour plots in the IRAC bands 
and in each case it is clear that the models are successful in 
matching the data, which again comes from Wang et al. 
(2010) and is cut to AT < 21. Figs. [5] and [6] indicate that by 
taking simple empirical spectra and evolving them back to 
higher z it is possible to get a reasonably accurate prediction 
for the average mid-IR colours of galaxies at a given redshift. 
Using colours to select galaxies in a specified redshift range 
has become a powerful tool, for example in the selection of 
Lyman-break galaxies (LBGs; e.g. Steidel et al., 2003; Adel- 
berger et al., 2004, Bielby et al., 2010), distant red galaxies 
(DRGs; Franx et al., 2003), extremely red objects (EROs; 
e.g. Roche et al., 2002), pBzK and sBzK galaxies (Daddi et 
al., 2004), and "BM/BX" galaxies (Steidel et al., 2004). Our 
PLE models may provide a simple underlying framework in 
which these colour selections can be better understood. 

The [3.6] - [4.5] : [5.8] - [8.0] colour-colour plot (Fig. 
[n];) has been used widely to select populations at different 
redshifts (Stern et al., 2005; Yun et al., 2008; Dey et al., 
2008; Devlin et al., 2009). The data in this plot visibly sepa- 
rate into two 'tiers' with the upper tier, redder in [3.6]-[4.5], 
hosting higher-z objects. Our models predict that the up- 
per and lower tiers should divide galaxies above and below 
z — 1.3; this is consistent with the finding of Devlin et al. 
(2009) that it cut galaxies at z « 1.2, especially given that 
they found « 15% contamination of higher-z sources in the 
lower tier and lower-z sources in the upper tier. Stern et al. 
(2005) identified an 'AGN wedge' in the upper tier, a region 
they suggest is populated by active galactic nuclei. For ref- 
erence we have marked this region in Fig. ^]p; it is located 
redward of the high-z end of our galaxy models in [5.8]-[8.0]. 
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Figure 7. Differential number counts in IRAC's 4 channels: (a) 3.6iJ.m, (b) 4.5|^ni, (c) 5.8(i.m and (d) S.On-m. The blue circles are data 
from Fazio et al. (2004) and the red squares are our own counts from public GOODS-N data. The solid black line shows the prediction of 
our model. The dot-dashed and dashed black lines indicate the contributions from late- and early-type galaxies, respectively. The green 
dotted line shows the combined contribution of these populations but from the stellar light only, i.e. without the emission from dust. 



Essentially, simple colour-redshift tracks present a vi- 
able way for observers to select galaxies in a given z-range. 
We find that sources which lie on the upper [z > 1.3) tier 
of the [3.6] - [4.5] : [5.8] - [8.0] plot also lie around the 
1.3 < z < 3 track in, for example, a. B — R : R ~ I oi 
R — I : I — K plot, so we emphasise that colour selection 
techniques in the context of PLE models or otherwise can be 
used successfully and consistently all the way through the 
optical, near-IR and out to mid-IR wavelengths. 

In summary, we have demonstrated that mid-IR colour- 
redshift data show no contradiction or inconsistency with 
the predictions of optically-defined galaxy evolution models 
with r = 1 Gyr for early-types and r — 9 Gyr for late-types. 



5 MID-IR NUMBER COUNTS 

5.1 IRAC bands 

Fig. [7] shows the number counts of galaxies at 3.6, 4.5, 5.8 
and 8.0[xm. The predictions of our model are compared to 
data from Fazio et al. (2004), who combine observations 
from three regions ranging from an ultradeep 5' x 10' field 



to the wide 3° x 3° Bootes field, in order to achieve good 
sampling of both the faint and bright end of the counts. We 
also compare to counts of GOODS-N sources by using SEx- 
tractor on the public Spitzer imaging of the that field (Dick- 
inson et al., in prep.), which agree well with Fazio et al.'s 
measurements. In both datasets stars have been excluded. 

Broadly speaking, the agreement between the model 
and the data is good. Although the model appears to over- 
predict the counts at the faint end in panels (a)-(c), the 
data here suffer from confusion. The fit at 5.8iJ.m is poorer 
than in the other bands, possibly due to our somewhat unso- 
phisticated modelling of the PAH features, which may have 
a more noticeable effect at 5.8iJ.m than elsewhere. While a 
more refined model of these features may improve the pre- 
diction, we do not consider this adaptation to be necessary 
for the purposes of this paper. Indeed, we note that the fit 
to the 5.8vim data is very similar to that achieved by the 
ACDM hierarchical formation models presented by Lacey et 
al. (2008). 

The plots also show the contribution made in each band 
by stellar light only (as opposed to dust). At 3.6vi.m and 
4.5iJ.m the stellar component accounts for the almost all of 
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Figure 8. The fraction of the galaxy number counts accounted 
for by dust emission, rather than stellar light, in each of the IRAC 
channels, based on the models in Fig. [7] The shift from starlight- 
dominated to dust-dominated emission appears to occur around 
A Ri 5— 6jim. 



the counts, while in the longer wavelength bands it yields 
an underprediction. In Fig. [8] we show the fraction of the 
total predicted counts arising from dust emission in each of 
the IRAC bands, clearly demonstrating the transition from 
starlight-dominated to dust-dominated galaxy counts as we 
move from the near- to the mid-infrared. 

In each band the very faintest source counts remain 
starlight-dominated, even at 8(i.m. The model suggests that 
direct stellar emission accounts for only ~ 10% of the 8jxm 
number counts at magnitudes brighter than tti = 13 (a flux 
of S'sum ^ 0.4 mjy), but makes up the dominant fraction 



18 {Si 



8|J.m 



< 4 



(> 90%) of Sfim sources fainter than m 

Overall, we have shown that a basic PLE model, orig- 
inally designed to match optical B-band counts, continues 
to provide a reasonably good description of the observed 
source counts out to the mid-IR wavelength regime at Sjim 
if a modest amount of dust absorption and re-emission for 
spiral galaxies is included. 



5.2 24 microns 

In Fig.[5]we show 24)i.m number counts, with data from Pa- 
povich et al. (2004) and our own source counts from public 
GOODS-N 24^m data (Chary et al., in prep.). We find that 
our model of normal galaxies, using the same dust param- 
eters which had reasonable success at 3.6 < A < 8.0 [im 
underpredicts the number counts by a factor of 4-5 (dashed 
line in Fig. [5|, accounting for only 22% of F24 > 30 viJy 
sources. This, then, is the wavelength regime where our ba- 
sic model first breaks down. 

However, we find that a relatively simple alteration to 
the model addresses this problem somewhat. We introduce 
a subset of our spiral galaxies, modelled on the dust features 



of Mrk33, as presented in the SINGS spectroscopic survey of 
local galaxies. This survey uses two-component dust models 
to match observed mid-IR SEDs of local galaxies (Draine et 
al., 2007). In most cases the dust emission is dominated by 
dust in the diffuse ISM, but some galaxies' dust emission is 
dominated by a warmer dust component, attributable pri- 
marily to photo-disassociative regions (PDRs), where the 
dust is exposed to a stronger intensity of radiation. These 
sources are found to have typically higher star-formation 
rates. 

The PDR dust is hotter and thus its emission peaks at 
shorter wavelength (see Fig. [2]). As a result of this strong 
PDR contribution, therefore, these galax;ies exhibit a signif- 
icant excess in flux in the mid-infrared. We flnd that they 
can make a signiflcant contribution in the A ^ 24|j.m regime, 
without adversely affecting the flt of the model at shorter 
optical to mid-IR wavelengths. 

18% of the 65 SINGS galaxies discussed in Draine et 
al. (2007) are of this type; these include Mrk33, NGC2798, 
NGC3049 and Tol89. We modelled an SED with such an 
infrared excess, based on the starburst galaxy Mrk33. Our 
SED is based on the observed mid-IR spectrun^Q of Mrk33 
from SINGS observations with IRS as well as broadband 
photometry from Dale et al. (2005) to constrain the dust 
blackbody at longer wavelength. This SED and the broad- 
band data appear in Fig. [S] We substituted this SED for a 
comparable fraction (16%) of the dusty late-type galaxies 
in our model by replacing the Scd galaxies (Table [TJ with 
this starburst population. This component is shown as the 
dash-dot line in Fig. [9] 

We find that including these sources in our PLE model 
yields a significant improvement in the fit to the counts, 
particularly at the bright end, suggesting that at 24fi.m the 
counts are dominated by sources with warmer dust and 
higher SFRs - primarily starburst galaxies. 

While the bright end of the number counts are ex- 
tremely well matched, at S2i\x-m < 0.5 mJy the counts are 
still underpredicted by a factor of a few. This suggests that 
our model is deficient in the number of predicted 24(i.m 
sources at high redshift, and comparing the predicted red- 
shift distribution of the model to the observed n{z) (Perez- 
Gonzalez et al., 2005) confirms that this is the case. 

Babbedge et al. (2006) suggest that while the galaxy 
luminosity functions at IRAC wavelengths show relatively 
little evolution, there is strong evolution in the 24jj-m LF, 
particularly out to z ^ 1. Including our Mrk33 component 
brings our model in line with local observations and in do- 
ing so enables us to match the bright-flux, low-redshift end 
of the counts. Allowing this warm-dust, high-SFR compo- 
nent to evolve more strongly could make it possible to ac- 
curately fit the 24|j.m counts across the full range of flux; 
this would be in line with the approach taken in previous 
analyses such as that of Babbedge et al. (2006) and Rowan- 
Robinson (2001; 2009). However, we explore here the alter- 
native possibility that the counts can be matched by ex- 
ploiting the negative 24jim fc-correction, combined with the 
moderate PLE which was successful at IRAC wavelengths, 
instead of invoking much stronger evolution. 

This is motivated initially by the finding that the shape 



^ available for download at |http://sings.stsci.edu7] 
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Figure 9. Differential number counts at 24|im. Blue circles 
are data from Papovich et al. (2004), green triangles are from 
Bethermin et al. (2010) and red squares are our own counts from 
public GOODS-N MIPS data. The lower solid line is the predic- 
tion of our model, with the 'normal' spiral galaxy contribution 
shown as a dashed line, Mrk33 galaxies shown by a dot-dashed 
line and the negligible early-type contribution shown as a dotted 
line. Mrk33 galaxies clearly dominate the counts, with the nor- 
mal galaxy component underpredicting the data by a factor of 
~ 5. The bright end is extremely well matched, but at the faint 
end the model is deficient by a factor of a few. We note that the 
spirals match the shape of the counts well, so we show also the 
prediction of the model where instead of Mrk33 we use an SED 
with the bright-end (z = 0) normalisation of Mrk33, but the k + e 
evolution of the normal spirals. This is the upper solid line. 



of the predicted source counts for the normal spirals matches 
weU the shape of the data, showing an upturn around 0.5-1.0 
mjy where the counts briefly increase at a super-Euclidean 
rate (Fig. |9]). The Mrk33 model, although matching the 
bright end normalisation well, has a much flatter shape with 
no increase at the faint end. 

The faint-end increase in the spiral model arises from 
a strong negative fc-correction, whereas in the Mrk33 SED 
the fc-correction is positive. As the SEDs are redshifted, the 
PAH region moves into the 24(im band. As this happens, the 
cool-dust spiral model becomes brighter (giving a negative 
fc-correction) but the Mrk33 model, which has a warm dust 
component, becomes fainter overall. (See Fig. [2]). 

Therefore, where we previously used the Mrk33 compo- 
nent, we now use instead a population with the same z = Q 
LF as the Mrk33 population (i.e. the same K — [24] colour), 
providing the good fit to the bright end counts, but with 
the k + e evolution of our normal spiral galaxies. We find 
this gives a good fit to the data across the full flux range 
(Fig. |9]). We are essentially using a hypothetical SED which 
has the same shape as the rest of our spiral galaxies but a 
redder K — [24] colour and a greater total IR luminosity. 

Using the appropriate bright end normalisation, then, 
PLE yields a good fit across 3 decades of 24)i.m flux. Since 
there is a decoupling of the stellar and dust emissions, this 



arbitrary normalisation is not unphysical, however unlike 
our initial inclusion of the Mrk33 component it is not obser- 
vationally motivated. 

We conclude that the 24(i.m band is the first regime 
where our simple bimodal PLE model begins to break down. 
Even with the Mrk33 galaxies included, the model does not 
have the sufficient ingredients to match the faint-flux, high- 
redshift counts. However, we have shown that these data can 
be well matched if one component of our spiral galaxies is 
altered to be much redder in A" — [24], although the data 
can alternatively be matched by aUowing for much stronger 
evolution, above the level found in our PLE model. 



6 FAR-IR/SUB-MM NUMBER COUNTS 

In Fig. [TO]we show number counts in 6 wavebands in the far- 
IR and sub- mm regimes. Counts at GO^m are from wide-field 
IRAS surveys (Lonsdale et al., 1990; Rowan- Robinson et al., 
1991; Gregorich et al., 1995), 70ti.m and 160iJ.m data come 
from more recent deep Spitzer surveys (Dole et al., 2004; 
Bethermin et al., 2010) and counts at 250-500vimcome from 
BLAST (Patanchon et al., 2009) and Herschel (Oliver et al., 
2010; Dunlop et al., 2010) observations. 



6.1 Galaxy model predictions 

In the far-IR (upper three panels) there appears to be a 
continuation of the result we found at 24(j.m: the galaxy 
model (solid black line), which continues to be dominated 
by the Mrk33 starburst component (dot-dashed black line), 
matches the bright end of the counts well, but perhaps be- 
gins to underpredict the faint end - this can be seen in the 
faintest 70vim counts. Our result at GOfim is consistent with 
that of Busswell & Shanks (2001), who also achieved a good 
flt to the 60[im counts using a similar model. 

While the flt to the 60-70ti.m counts is, overall, rea- 
sonable, the model does very poorly at matching IGOvim 
counts, with even the bright end underpredicted. Moving 
to submillimetre wavelengths (lower three panels) , although 
the bright end appears to match well, at faint fluxes the 
model underpredicts the data by as much as 1-2 orders of 
magnitude. Therefore, normal galaxies which dominate ob- 
servations at optical to mid-infrared wavelengths account for 
only a few percent of sources seen in the sub-mm. 



6.2 Obscured AGN model predictions 

The conventional view of the submillimetre-bright galaxy 
population observed at 850[Am is that they are high-z 
ULIRGs with very high levels of star-formation fuelling their 
extreme bolometric luminosities (Small et al., 1997; Barger 
et al, 1998; Alexander et al., 2005). However the possibil- 
ity remains that obscured AGN make a major contribution 
to this population (Almaini et al., 1999; Gunn & Shanks, 
1999): in a recent paper (Hill & Shanks, 2010) we demon- 
strated that a model of obscured AGN could give a good 
flt to the observed 850vim source counts and extragalactic 
background. We showed that the predictions of this model 
for the sub-mm contribution made by X-ray-detected AGN 
were in line with our measurements for X-ray sources in the 
Extended Chandra Deep Field South. 
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Figure 10. Differential number counts at 60—500 microns. We show the prediction of our galaxy model (solid black line) as well as its 
constituent components - normal spirals (dashed black) and the Mrk33 component (dot-dashed black). This model gives a reasonable 
fit to the bright-end data (except at 160|^m) but underpredicts the faint end by an increasingly large factor. We show also the predicted 
contribution from an obscured AGN model (long-dashed red) and the combined galaxy-f AGN prediction (solid red), which dramatically 
improves the fit at long wavelengths. 
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Figure 11. Redshift distributions at 250n.m, 350|J.m and 500|J.m, showing our model prediction (for galaxies and obscured AGN combined; 
solid line) against the observed distributions published by Chapin et al. (2010; dashed line). The histograms have been normalised to 
the same total number of sources. In all three bands there is good agreement between the model and the data. 



The model we use is that of Gunn & Shanks (1999), 
consisting of a obscured AGN enshrouded with a fiat distri- 
bution of column densities between A'h ~ 10^®'^ cm~^ and 



10^ 



cm and dust masses inferred from a Galactic gas to 



dust ratio. This non-unified AGN model is known to match 
measurements of the X-ray background (Gunn, 1999). As 
with our galaxy model we assume a dust temperature of 
Td = 30K for the AGN and all radiation absorbed by the 
obscuring medium is assumed to be isotropically re-radiated 
in the infrared. The population evolves through PLE only. 

This model gives a good fit to 850jtm number counts 



(Hill & Shanks, 2010). Although it underpredicts the faint 
end, at this point our model of normal galaxies makes up 
the counts well. We now apply the obscured AGN model to 
shorter wavelengths, keeping all parameters the same as in 
our previous paper. In Fig. 1101 the contribution from the 
AGN is shown as a red long-dashed line, and the combined 
AGN and galaxy model is shown as a solid red line. At 60(xm 
and 70iJ.m, the AGN are entirely negligible across the full 
flux range. At IGOjtm, their contribution is roughly equal 
to that of the galaxies, however the combined galaxy-AGN 
model still falls short of the number counts. 
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At 250-500 [im, we find that the AGN model can dra- 
matically improve the fit to the data. At 250|i.m and 350n.m, 
the prediction at faint fluxes (S < 30 mjy) is still lower 
than the data, but sources here are beyond the confusion 
limit (Nguyen et al., 2010) and so the data must be treated 
somewhat cautiously. At 500[im, the bright end is a little 
overpredicted by the AGN model, however we suggest that 
overall the model matches the sub-mm data with reasonable 
success, and we emphasize that no aspect of the model has 
been altered to fit these source counts. 

In Fig. Illlwe show predicted sub-mm redshift distribu- 
tions compared to data from a BLAST survey (Chapin et 
al., 2010). The agreement between the model and the data 
is generally very good. Dunlop et al. (2009) show that the 
250nm N{z) is bimodal, with the population divided broadly 
into low-z spirals and high-z submillimetre galaxies. This is 
replicated in our model - the low-redshift peak arises from 
our spiral galaxies and the z > 0.5 sources are made up of 
the obscured AGN. 

In summary, our PLE galaxy model continues to per- 
form reasonably well out to A ~ 70(i-m, but is unable to 
match source counts at longer wavelengths. At this point we 
invoke the PLE model of obscured AGN which we have pre- 
viously shown to match 850 nm data, and find that it also 
gives a reasonable fit at 250-500 (am. 



7 DISCUSSION 

We have shown that our basic PLE galaxy model can give 
a good description of galaxy populations out to 8(i.m, with 
starlight dominating the counts out to « 6vim and dust emis- 
sion dominating after that out to 8 microns. With a star- 
burst component included in line with local observations, the 
model can continue to match observations out to « 70 [im 
with reasonable success. A caveat is that the 24-70 micron 
counts have been fitted here with a K-correction that comes 
from a hybrid galaxy template, formed from a stellar black- 
body and a PAH component in a ratio not found in local 
galaxies; others have interpreted data at these wavelengths 
as requiring strong evolution and this remains a viable al- 
ternative interpretation. Either way, the success of the PLE 
model in fitting data at these long wavelengths is still im- 
pressive, given that PLE represents the simplest approxima- 
tion other than no-evolution models. 

There is an increasing consensus surrounding the hier- 
archical galaxy evolution scenario, in which galaxies form 
at late times through a gradual build-up of stellar mass via 
the merging of smaller systems. This motivates the ques- 
tion as to whether PLE models should be treated simply as 
a convenient phenomenological description of galaxy popu- 
lations, which can match - and, more valuably, predict - 
observations, or whether they may in fact have a real, phys- 
ical basis. In the latter case, the implication would be that 
galaxies form at high-z in a monolithic collapse event, and 
passively evolve thereafter. But if these models represent a 
phenomenological description only, then it must be that hi- 
erarchical galaxy evolution is able to present the appearance 
of generally passive evolution. 

The detection of massive galaxy populations at z > 2 
has proven a challenge for hierarchical models (Baugh et al., 
1998; Somerville, Primack & Faber, 2001; Bower et al., 2006) 



since hierarchical formation should intuitively produce few 
massive structures at high-redshift given that the process of 
formation is one of gradual build-up. However, it has been 
shown that massive galaxies at z > 2 can be explained by 
invoking AGN feedback (e.g. Bower et al., 2006) and/or cold 
gas accretion (Dekel et al., 2009; Kang et al., 2010). At lower 
redshift, z < 1, hierarchical models predict a rapid evolution 
in the number of massive systems, with perhaps a tenfold 
drop in the number of massive {Ad > lO^^M©) galaxies be- 
tween z = 1 and the present (Baugh et al., 2003). Observa- 
tions, on the other hand, have shown no evidence for density 
evolution of massive ~ L* ellipticals out to z = 0.9 (Wake et 
al., 2006; Brown et al., 2007; Cool et al., 2008). However, it 
has been possible to reconcile the slow evolution of observed 
early-type galaxies with the hierarchical model by use of a 
suitable Halo Occupation Distribution (HOD). 

PLE models are entirely consistent with the presence of 
massive galaxies at high-z and the small amount of evolution 
z < 1, however these models have had only mixed success 
in matching the clustering of massive galaxies. Wake et al. 
(2008) suggest that the evolution in the clustering of lumi- 
nous red galaxies (LRGs) between z = 0.6 and z = 0.2 rules 
out passive evolution, particularly at small scales, finding 
that a hierarchical model with a low merger rate provides 
a better fit. However, Sawangwit et al. (2009) show that a 
passive model, with no change in comoving source density, 
is fully compatible with their clustering measurements in 
a much larger LRG sample. Thus PLE has not only been 
found to be an excellent phenomenological fit to the data, it 
is still possible that it may even represent a good physical 
model as well. 

At wavelengths of 160vim and longer our galaxy PLE 
model finally fails. We then invoked the non-unified, ob- 
scured AGN model of Gunn & Shanks (1999) which fits 
the 870 micron source counts at bright fluxes, leaving the 
galaxy model to contribute at fainter fluxes (Busswell & 
Shanks, 2001; HiU & Shanks, 2010). We showed that this 
model also fits the Herschel counts at 250, 350 and 500iJ.m. 
Given that this could remove the need for high luminosity, 
massive starbursts to fit the sub-mm counts, then this AGN 
model could be said to ease the problem that such massive 
galaxies at high redshifts represent for the standard ACDM 
cosmological model. There would be for example no need 
to invoke a top heavy IMF to explain such sources (Baugh 
et al., 2005). It would only remain to check whether the 
model could accommodate the galaxy PLE phenomenology 
at shorter wavelengths (A < 160vim). 

The AGN model we use is non-unified, with absorbed 
AGN having brighter sub-mm fluxes than unabsorbcd AGN 
(see also Page et al., 2004). Hill & Shanks (2010) suggested 
that if a unified rather than non-unified AGN model was 
applied, the AGN contribution to the sub-mm would drop 
dramatically and another component would be required to 
fit the sub-mm source counts. Conventionally, this is taken 
to be a population of highly evolved starburst galaxies with 
the bolometric luminosities of QSOs. We now consider this 
alternative possibility in the context of a PLE model where 
these sources are modelled as dust-obscured early-type pro- 
genitors. We have therefore modelled a population of early- 
types that evolve according to a r = 1 Gyr model with a 
Chabrier IMF and dust absorption in the range Ab = 0.3-1 
mag, which re-radiates optical light into the FIR assuming 
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the usual 30K dust temperature. Although this model pro- 
duces a K band evolution of 2-3 mag at 2 « 2-6 relative 
to the a; = 3 IMF model or the fc-correction, we found the 
contribution to the sub-mm counts to be negligible. 

In addition, the Chabrier IMF produces strong evolu- 
tion at 2 ~ 1 and would need a dust absorption of y4s ~ 2 
mag to avoid over-predicting the high-« tail of the A" < 20 
redshift distribution. This level of extinction would destroy 
the viability of the PLE models to fit galaxy colour and 
count distributions through the optical and NIR bands. Ap- 
pealing to a top-heavy IMF might still enable the sub-mm 
counts to be fitted, but such an IMF could not be allowed to 
apply at 2; « 1-2 because the problem of overpredicting the 
K20 n{z) would be exacerbated. Somewhat paradoxically, 
although the semi-analytic model that invokes a top-heavy 
IMF (see Baugh et al., 2005; Lacey et al., 2008) fits weU 
the sub-mm counts, it appears to seriously underpredict the 
average K20 redshift (Gonzalez-Perez et al., 2009). These 
authors also show that the model of Bower et al. appears to 
fit the K n{z) better, but fails to fit the sub-mm counts. 

Certainly in the context of our PLE model there is 
greater observational motivation for a dwarf-dominated IMF 
than a top-heavy IMF: the top-heavy IMF is not required 
since sub-mm counts are matched by the AGN population, 
and the dwarf-dominated IMF can therefore provide an ex- 
planation for the absence of a detectable initial burst of star- 
formation from early-type galaxies at very high redshift - 
this burst may be absent due to a lack of bright stars in the 
IMF of early- types. Arguments against such a model include 
the apparently prompt enrichment of early type galaxies 
evidenced by their a-element enhancement (Tinsley, 1979; 
Trager et al., 2000). But the excellence of the fits of our sim- 
ple galaxy count models from the UV through the K band to 
the sub-mm means that this alternative idea for the origin 
of early-types is still worthy of serious consideration. 



underprediction of the faint-end counts. If we instead use 
a population of spirals which have our standard late-type 
SED but the z = K - [24] colour of Mrk33 (essentially, 
the normal spiral SED normalised to brighter fiux), we fit 
the counts across the full fiux range. 

• At sub-mm wavelengths (A = 250, 350, 500vim) the 
galaxy PLE model fails - we find that normal galaxies ac- 
count for only a few percent of sources here. However, we 
show that a PLE model of obscured AGN, which we have 
previously shown to give a good fit to observations at 850(i.m 
(Hill & Shanks, 2010), can again give a good fit to the counts 
and redshift distributions in these bands. 
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